The Mass and Leptonic Decay Constant of 
Z} s o(2317) Meson in the framework of thermal 

QCD sum rules 

2 ! El§en Veli Veliev *, Giil§ah Kaya ** 

Physics Department, Kocaeli University, Umuttepe Yerle§kesi 

a '■ 41380 Izmit, Turkey 

■ * e-mail: elsen@kocaeli.edu.tr 

CN ' ** e-mail: gulsahbozkir@kocaeli. edu.tr 

(N : 

i-C ■ Abstract 

In the present work, we assume D s o(2317) meson as the cs state and study its param- 
eters at finite temperature using QCD sum rules. It is calculated the annihilation and 
r~| | scattering parts of spectral function in the lowest order of perturbation theory. Taking 

into account perturbative two-loop order a s corrections and nonperturbative corrections 
— . up to the dimension six condensates it is investigated the temperature dependences of 

mass and leptonic decay constant of D s q(2317) meson. 

t> ■ 

1 Introduction 

In 2003 BaBar Collaboration discovered a positive-parity scalar charm strange meson D s o(2317) 
with a very narrow width fl], which was confirmed by CLEO Collaboration [2] and BELLE 
Collaboration |3J later. This observed state has attracted much attention because its measured 
mass and width do not match the predictions from potential-based quark models [4j. To resolve 
the difficulties, many theoretical models have appeared in the literature. Various theoretical 
models, based on the cs quark structure, are suggested to explain the low mass and the narrow 
width for the D s0 (2317) meson [5]-[10]. QCD sum rule analysis in [TT] . [T2] supports the cs 
postulation of nature _D s o(2317). Apart from the quark- ant iquark interpretation, this state has 
been interpreted as a DK molecule [13], a D S 7T molecule [14] . a csqq four-quark state [15], and a 
mixing of the conventional state and the four-quark state [16]. Also this state was investigated 
in the framework of chiral symmetry considerations [TTj . 

In this work, we assume -D s o(2317) meson as the cs state and study its parameters at finite 
temperature using QCD sum rules [18]. The extending of QCD sum rules method to finite tem- 
perature has been made in the paper [19] . This extension based on two basic assumptions, that 
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the Operator Product Expansion (OPE) and notion of quark-hadron duality remain valid at 
finite temperature, but the vacuum condensates must be replaced by their thermal expectation 
values. The thermal QCD sum rule has been extensively used for studying thermal properties 
of both light and heavy mesons as a reliable and well-established method [20]-|24j . 

In the present work, we calculated the temperature behavior of mass and leptonic decay 
constant of D s0 (2317) meson. The knowledge of leptonic decay constants is needed to predict 
numerous heavy flavor electroweak transitions and to determine Standard Model parameters 
from the experimental data. Also leptonic decay constants play essential role in the analysis of 
CKM matrix, CP violation and the mixings BjBd, B S B S . 

This paper is organized as follows. In section 2, we calculated the annihilation and scattering 
parts of spectral density and give the expression for the perturbative scalar spectral function 
up to two-loop order a s corrections. Also nonperturbative contributions up to the dimension 
six condensates [16] are summarized. Section 3 contains our numerical analysis of the mass and 
leptonic decay constant using Borel transform sum rules. 

2 Thermal QCD sum rule for the scalar charm strange 
meson 

The starting point for the sum rule analysis is the two-point thermal correlator 



where J(x) = (m c — m s ) : s(x)c(x) : is heavy-light quark current and has the quantum numbers 
of the -D s o(2317) meson, m c and m s are charm and strange quark masses respectively, s quark 
mass is not neglected throughout this work. Thermal average of any operator O is determined 
by following way 



where H is the QCD Hamiltonian, and /3 = 1/T stands for the inverse of the temperature T 
and traces carry out over any complete set of states. According to the basic idea of the QCD 
sum rule, we must calculate this correlator in terms of the physical particles (hadrons) and in 
quark-gluon language, and then equate both representations. First let us calculate theoretical 
part of the correlator Eq. (1). Up to a subtraction polynomial, which depends on the large q 2 
behavior, II(g 2 ) satisfies following dispersion relation [19], [21], [22] 
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where p(q) = ~ImH(q) tanh(^ a ) is spectral density. In order to calculate two-point thermal 
correlator in the lowest order of perturbation theory we use quark propagator at finite temper- 
ature [25] 

S n {q) = (7% + m)(- \— + 2mn(u q )5(q 2 - m 2 )), (4) 

K q z — m z + is 



Here n(u q ) is the Fermi distribution function, n(u q ) = [exp(/3oj q ) + 1] 1 and u q = y/q 2 + m 2 . 
After some calculations we find that perturbative part of spectral density is given by 

< t\ f rfk ^i - k 2 + k • q - LOiq + rn c m s 

Ppert{q,-L > — / ,„ no 

J [ZTiy U)\UJ 2 

x [(1 - m - n 2 )5(q - ui - u 2 ) + (ni - n 2 )S(q - wi + w 2 )]. (5) 



Here cui = ^cf + m 2 . and co>2 = y (k — q) 2 + m 2 . . Note that spectral density involves two 
pieces, one is called the annihilation term, p a , P ert{s,T), which survives at T = 0. Other term 
is called scattering term, p s>per t(s,T), which vanishes at T = 0. As can be seen, delta function 
5(qo — u>i — uj 2 ) in Eq. (5) gives the first branch cut, q 2 > (m c + m s ) 2 , which coincides with zero 
temperature cut and describes the standard threshold for particle decays. On the other hand, 
delta function 5(qo — uo\ + u 2 ) in Eq.(5) shows that an additional branch cut arises at finite 
temperature, q 2 < (m c — m s ) 2 , and this new branch cut corresponds to particle absorption from 
the medium. Therefore, delta functions 6(q —uii —u) 2 ) and S(q —U\ +uj 2 ) in Eq.(5) contribute 
in regions (m c + m s ) 2 + q 2 < q$ < 00 and < q 2 < q 2 + (m c — m s ) 2 respectively. Taking into 
account these contributions the annihilation and scattering parts of spectral density in the case 
q = can be written as 

p a , P ert{s,T) = p {s)[l - n(^-(l + — s -)) -n(-^-(l ( 6 ) 

Ps, per t(s,T) = p Q ( s )[n[y^(l + 6 - s )) -n{- ^-(l - c s s ))\, (7) 

Here po{s) is the correlation function in the lowest order of perturbation theory at zero tem- 
perature and given by 

Po(s) = 3{mc - ms)2 q\s)v%s), (8) 

where q(s) = s — (m c — m s ) 2 and v(s) = (1 — Am s m c /q(s)) 1 ^ 2 . The contribution of perturbative 
two-loop order a s corrections to the spectral density in perturbation theory at zero temperature 
can be written as 1261: 



An 

Pi(s) = ^Po(s)f(x), (9) 



where x = m 2 /s, a s = a s {m 2 ) and 



9 3 x 

fix) = - + 2Lz 2 (a;) + lna;ln(l-x) \n(l/x- 1) -ln(l -x) +x\n(l/x- 1) \nx. (10) 

4 2 1 — x 



Here Li 2 (x) = — f£ dt Xa ^ l ~^ is dilogarithm function. Note that in a s corrections terms the 
strange quark mass is set zero, though in numerical analysis, the mass of the strange quark is 
taken account. The subtraction terms in Eq. (3) are removed by using the Borel transformation, 
therefore we will omit these terms. The non-perturbative contribution at zero temperature to 
the correlator has following form 
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m c A(0|ss|0) 



1 - ^r(3 - A) - \e 2 {\ - A) + ^ 3 (1 + A - 4A 2 + 2A 3 ; 



+ 



12tt 
M 2 



A(0|a s G 2 |0) 



3e 1- -A 
3 



2A 2 -2A(1-A) ln(eA) 



f^(0|ss|0)A 2 (l - A)(l - e(2 - A)) - 4-^a s (0|ss|0) 2 A 2 (2 - A - A 2 ) 



2m. 
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which arises in the framework of the OPE and parameterized by vacuum expectation values of 
quark and gluon fields in the QCD Lagrangian. In Eq. (11) A = m 2 c /(Q 2 + m 2 ), e = m s /m c 
and terms are organized according to their dimension. It is assumed, that the expansion (11) 
also remains valid, but the vacuum condensates must be replaced by their thermal expectation 
values [15] . For the light quark condensate at finite temperature we use the results of [27j,[28j 
obtained in chiral perturbation theory and temperature dependence of quark condensate in a 
good approximation can be written as 



<0|gg|0>[l - 0.4(— ) - 0.6( x 



;i2) 



where T c is critical temperature. The low temperature expansion of a gluon condensate is 
proportional to the trace of the energy momentum tensor [25] and can be approximated [23] as 



(a s G 2 ) = (0|a s G 2 |0) 
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(13) 



Also, we have used for the mixed condensate the parameterization 

g(qa^G^q) = M 2 (qq) 



(14) 



and deduced the value of the QCD scale A from the value of a s (M z ) = 0.1176 . 

Our next task is the calculation of the physical part of the correlator (1). According to the 
basic idea of quark-hadron duality assumption, the right-hand side of eq. (1) can be evaluated 
in a hadron-based picture. Equating OPE and hadron representations of correlation function 
and using quark-hadron duality the central equation of our sum-rule analysis takes the form: 

f 2 {T)m\T) _ pCT) 

'(m c +m s ) 2 



Q 2 + m 2 (T) 



dc Pa,pert{s,T) + pl{s) 



^ds^^+u^n (is) 



s + Q 2 
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where / and m are the leptonic decay constant and mass of D s0 (2317) meson respectively. 
Note that in Eq.(15) the mass and leptonic decay constant were replaced by their temperature 
dependent values. The continuum threshold also depends on temperature; to a very good 
approximation its scales universally as the quark condensate 



so CO = s - 



(qq) 



+ (m c + m s 



<0|g?|0) v s 

where in the right hand side Sq is hadronic threshold at zero temperature: Sq = s(T = 0) 



(16) 



3 Numerical analysis of mass and leptonic decay con- 
stant 



In this section we present our results for the temperature dependence of D s0 (2317) meson 
mass and leptonic decay constant. Performing Borel transformation with respect to from 
both sides of equation (15) and taking the derivative with respect to 1/M 2 from both sides of 
obtained expression, and making some transformations we have 



m 2 (T) = B(T)/A(T), 



(17) 



where 
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A(T) = I ds(p atPert {s) + pi{s)) exp ( - 

J (m c +m s ) 2 v ' ^ 

f(rn c -m s ) 2 / S \ 

i dsp StPert (s) exp ( - —J + n np (M , T) 
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(19) 



n np (M 2 , T) = ml(ss)e-e [l - |e + Ue - /fe 2 (l - U) + \e\l + p - 2/3 2 + \fr 



1 /a s G 2 



m 2 e 



1 - 3e(l - -p + (3 2 - 2p(hi(pe) + 7 - 1) + (3 2 (HM + 
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12 \ 7T 

+^M 2 m c /3He^[l - i/9 - 2e(l - - ~/3 
-^7rpa s (s S ) 2 /3e^(12 - 3/3 - /3 2 ), (20) 

ol 

where /3 = m 2 /M 2 and B(T) = -m 2 c ^p. 

For the numerical evolution of the above sum rule, we use QCD impute parameters showed 
in Table 1. The criterion we adopt here is to fix in such a way so as to reproduce the zero 
temperature values of meson mass and leptonic decay constant. D s q(2317) meson mass as a 
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Table 1: QCD input parameters used in the analysis. 



Parameters References 




function of temperature are shown in Fig.l, Fig. 3 and Fig. 5 at continuum threshold values 
so = 7.5; 8.0; 8.5 GeV 2 respectively. As seen, mass decreases with increasing temperature and 
mesons lose approximately 10—15 percent of its mass at T = 150 MeV temperature. The results 
for leptonic decay constants are shown in Fig. 2, Fig. 4 and Fig. 6 at continuum threshold values 
s = 7.5; 8.0; 8.5 GeV 2 respectively. As can be seen / decreases with increasing temperature and 
vanishes approximately at critical temperature. This situation may be interpreted as a signal 
for deconfinement and agrees with light and heavy mesons investigations [16], [23]. Numerical 
analysis shows that the temperature dependence of / is the same, when M 2 changes between 
1.5 GeV 2 and 3 GeV 2 at fixed values of continuum threshold. Obtained results can be used for 
interpretation heavy ion collision experiments. 
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Figure 1: Temperature dependence of meson mass at sq = 7.5 GeV 2 . 
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Figure 2: Temperature dependence of leptonic decay constants at sq = 7.5 GeV 2 . 
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Figure 3: Temperature dependence of meson mass at sq = 8.0 GeV 2 . 
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Figure 4: Temperature dependence of leptonic decay constants at sq = 8.0 GeV 2 . 
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Figure 5: Temperature dependence of meson mass at sq = 8.5 GeV 2 . 
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Figure 6: Temperature dependence of leptonic decay constants at sq = 8.5 GeV 2 . 
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